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We introduce an operational procedure to determine, with arbitrary probability and accuracy, 
optimal entanglement witnesses for every multipartite entangled state. This method provides an 
operational criterion for separability which is asymptotically necessary and sufficient. Our results 
are also generalized to detect all different types of multipartite entanglement. 
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Entanglement, first noticed by Einstein, Podolsky, and 
Rosen [1], is at the heart of quantum mechanics. Quan- 
tum teleportation, superdense coding and cryptogra- 
phy [2] are achieved only when one deals with insep- 
arable states. Thus, the determination and quantifi- 
cation of entanglement in a composite quantum state 
is one of the most important tasks of quantum infor- 
mation theory. A finite-dimensional density operator 
Pi...n G B{Hi (g) ... (g) Hn) (the Hilbert space of bounded 
operators acting on i7i (g) ... (g Hn) is separable iff it can 
be written as a convex sum of separable pure states: 



Pl...n = ^Pi\'^i)ll{'4'i\ •■• (g \tpi)nn{i>i 



(1) 



where {pi} is a probability distribution and \^i)k are vec- 
tors belonging to Hilbert spaces Hk- Despite the sim- 
plicity of this definition, no operational necessary and 
sufficient criterion have been found for the separability 
problem until now. Moreover, it was shown by Gurvits 
[3] that this problem is NP-HARD. In this letter, we 
present a procedure to determine, with a chosen proba- 
bility, if a given state is entangled. In order to do that, 
we apply a class of convex optimization problems known 
as robust semidefinite programs (RSDP) to the concept 
of entanglement witness (EW) which we briefly recall. 

An operator pi...„ is entangled iff there exists a self- 
adjoint operator W € B{Hi g) ... g) Hn) which detects 
its entanglement [4], i.e., such that Tr{Wpi,,,n) < and 
Tr(VF (Ti...„) > for all (tab separable. This condition 
follows from the fact that the set of separable states is 
convex and closed in B{Hi (g ... g) Therefore, as a 

conclusion of the Hahn-Banach theorem, for all entan- 
gled states there is a linear functional which separates 
them from this set. We will deal in this paper only with 
normalized entanglement witnesses such that tr{W) = 1. 



Definition 1 A hermitian operator Wopt G B{Hi g 
Hn) is an optimal EW for the density operator pi.. 

Tr{WoptPi...n) < Tr{Wpi...n) 
for every EW W . 



(2) 



Although this definition of OEW is different from the one 
introduced in [5], the optimal EWs of both criteria are 
equal. 

We may now express the search of an optimal EW for 
an arbitrary state pi...„ in terms of a robust semidefi- 
nite program (RSDP) . A semidefinite program (SDP) 
consists of minimizing a linear objective under a linear 
matrix inequality (LMI) constraint, precisely, 

minimize c^x subject to 



F(x) = Fo 



x,F,, > 



(3) 
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where c e C" and the hermitian matrices Fi — F- e C"^" 
are given and x € C™ is the vector of optimization vari- 
ables. F{x) > means F{x) is hermitian and positive 
semidefinite. SDPs are global convex optimization pro- 
grams and can be solved in polynomial time with interior- 
point algorithms [6]. For instance, if there are m opti- 
mization variables and F{x) is a nxn matrix, the number 
of operations scales with problem size as 0{m?n^). SDPs 
have already been used in different problems of quantum 
information theory [7] and also in the separability prob- 
lem [8]. An important generalization of (3) is when the 
data matrices Fi are not constant, i.e., they depend of a 
parameter which varies within a certain subspace. This 
family of problems, known as robust semidefinite pro- 
grams, is given by: 

minimize c^x subject to 



F(x,A) = Fo(A)-f^a;,F,(A)>0, VAeP (4) 

i=l 

where !> is a given vectorial (sub)space. Note that prob- 
lem (4) is more difficult to solve than (3) , since one must 
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find an optimization vector x such that F{x, A) is pos- 
itive semidefinite for all A G P. One often encounters 
SDPs in which the variables arc matrices and in which 
the inequality depends affinely on those matrices. These 
problems can be readily put in the form (3) by introduc- 
ing a base of hermitian matric;cs for each matrix variable. 
However, since most of optimization solvers [9] admit 
declaration of problems in this most general form, it is 
not necessary to write out the LMI explicitly as (3), but 
instead make clear which matrices are variables. Equal- 
ity constraints involving the optimization variables can 
also appear in (3) and (4) without any further computa- 
tional effort. We can now enunciate the main result of 
this letter. 

Theorem 1 A state G B{Hi (g) ... ® H^) is entan- 
gled, i.e., can not be decomposed as (1), iff the optimal 
value of the following RSDP is negative: 

minimize Tr{Wpi,,,n) subject to 

dn dji dn dji 

(<-<-i«ii-«in-i (5) 

Wii...i„_iji...j„_J > 

Tr{W) = 1, Voi^ eC, l<fc<n-l 
where dn is the dimension of Hn, Wi^,,,i^_.^j^,,,j^_^ = 

l{i\ ^ ...^n-l {i\W\j}n-l ® ... ® \j)l e B{Hi ... Hn-l) 

and \i)k is an orthonormal base ofHk- If pi...n is entan- 
gled, the solution matrix W which minimizes Tr{Wpi,,,n) 
is the OEW for p\,,,n- 

proof: First wc have to show that (5) is a genuine 
RSDP. Note that Wi^,,,i„_iji...j„_i and the objective 
Tr{W pi,,,n) are both linear in the matrix variable W. 
Thus (5) can be put in the form (4), where V, in this 
case, is C''". A state p\...n is entangled iff there exists an 
operator W such that Tr{W pi,,,n) < and i{tp\ (g) ... (g)„ 
(V'|W|'0)„(g)...(g)|V')i > for aU states \tp)k G Hk- There- 
fore, the matrix i(i/'|(8'...(8'n-i(V'|W^|V')n-i<^---<^l^)i > 
has to be semidefinite positive for all \tp)k G Hk ■ Let- 
ting \tp)k = J2j (^j\j)k, where \j)k is an orthonormal base 
of Hk, it is straightforward to show that the optimal W 
given by (5) is the OEW of /9i...„. QED. 

In spite of the similarity between (3) and (4), RSDPs 
are in general very hard optimization problems. Actually, 
it was proved that robust semidefinite programs in the 
form of (5) are NP-HARD [10]. 

Corollary 1 The determination of the OEW for an ar- 
bitrary state pi...n is a NP-HARD problem. 

Since (5) is computationally intractable, it is natural to 

search for approximations of it in terms of SDPs, which 
are very efficiently solved. These relaxations of RSDP 



have been intensively studied [11] in the past years and 
can be classified as deterministic or probabilistic. In this 
letter we will focus on the latter, where one seeks a fea- 
sible solution to most of the possible values of the vary- 
ing parameters. The results of applying deterministic re- 
laxations to (4), which yields new separability sufficient 
criteria, was reported in [12]. Our methodology will be 
based on the concept of e-level solution introduced in [13] . 

Consider the most general form of RSDP given by (4) . 
Assume that the support V for A is endowed with a a- 
algebra and that a probability measure P over this alge- 
bra is also assigned. Let x G C™ be a candidate solution 
to (4). The probability of violation of x is defined as: 
V{x) = P{A e V : F{x,A) < 0}. For example, in (5), 
where the varying parameters are uniformly distributed 
over C''", V{x) measures the percentage of parameters 
such that the linear matrix inequality is violated. 

Definition 2 Let e G [0,1]. We say that a hermitian 
operator W is an e-level entanglement witness, e-W, if 

V{W) = P{aeS: Tr{Wa) < 0} < e 

where S is the subspace of separable density operators. 

The concept of optimal e-level entanglement witness is 
totally analogous to the one of definition (1), but now 
(2) has to hold for every e-level EW. The importance of 
this new class of hermitian operators is that, in contrast 
to the case of genuine EW, e-level optimal EW can be de- 
termined with a priori chosen probability in polynomial 
time for every multipartite state. 

Theorem 2 Let e G [0, 1], /? G [0, 1] and N > ^^^^ - 
1, where D is the dimension ofHi®...®Hn- Assume that 
N independent identically uniformly distributed samples 
aj^ , a^^ , aj^, 1 < ji < dn and 1 < I < n — 1, are drawn. 
Then the optimal e-EW for a state pi...n is given with 
probability at least 1 — (3 by the solution of the following 
semidefinite program: 

minimize TriyV p\,,,n) subject to 

dn dn dn dn 

EE- E E ((4r-«-.r4-<-. 

»l = ljl = l »n-l=l jn-l=l 

(6) 

Wi^...in-dl-3n-l) > 

Tr{W) = 1, 1 < fc < AT 
where dn is the dimension of Hn, — 

x{i\®...®n-l{i\W\j)n-l®.-®\j)l e B{Hx®...®Hn-l) 

and \ j)k is an orthonormal base ofTLk. 

proof: According to [14] , an e-level solution of a RSDP 
can be obtained with probability 1 — /3 from a sampled 

convex program, where the robust linear matrix inequal- 
ity is replaced by A^ > — 1 independent identically 
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distributed samples chosen according to probability P, 
where r is the number of optimization variables of the 
problem. The result follows in a straightforward manner 
if one notices that problem (5) has D{D + 1) optimization 
variables (the number of distinct real entries of W) and 
that P in this case is uniform. QED. 

Notice that theorem (2) gives a sufficient condition for 
separability, which is asymptotically also necessary. In 
fact, it is possible to determine, with any desired pre- 
cision and probability, if any state is entangled or not. 
Nevertheless, one must always consider the trade-off be- 
tween the accuracy of the results and computation effort. 
Although a priori feasibility levels are given by the for- 
mer theorem, the optimization problem yields in general 
much better results. Once a solution has been deter- 
mined, it is possible to make an improved estimate of the 
level of feasibility using Monte-Carlo techniques. In order 
to do that, generate a new set of N independent identi- 
cally uniformly distributed samples ttj^ and construct the 
empirical probability of constraint violation, Vemp{W) = 



Er=i ® ■■■ ®n {MW\iJj)n ® ••• ® |V^)i < 0), where 



N Z—'i=l 

is the indicator function. Then, the classical Cher- 
noff inequality guarantees that — VeTOp(W)| < e 

holds with confidence grater than 1 — /? , provided that 



TV > 



2e2 



(7) 



samples are drawn. Another important performance pa- 
rameter is the minimum eigenvalue over the violated con- 
straint. It can also be obtained empirically and it is very 
useful to determine if the solution obtained is accurate. 

We present now the first example for which we applied 
our techniques to determine an approximate optimal en- 
tanglement witness. We used MATLAB and the package 
SEDUMI [9] to implement and solve the SDP. Consider 
the Horodecki 3X3 bound entangled states [15]: 
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(8) 



where a S [0,1]. This family of states is particularly 
interesting because the Peres-Horodecki criterion fails to 
detect its entanglement, i.e. they are positive partial 
transpose entangled states. Using N = 1200 samples in 
each test, we were able to detect entanglement for all 
values of a, except for and 1. The expectation value of 
the OEW for each p(a) is shown in Fig. 1. The empiri- 
cal probability (Vernp(W)) of violation and the minimum 
eigenvalue over the violated constraint {Xmin), calculated 
using N = 10^ samples, were both negligible, showing 
that the algorithm converged. 



0.2 0.4 



FIG. 1: Tr{Woptp{a)) X a, for the 3 X 3 Horodecki bound 
entangled states. 



We have applied our methodology to a large number of 

2x2 and 2x3 states, namely. 5000 (five thousand) ran- 
dom states of each kind. Since in this case the positive 
partial transpose criterion [16] gives sufficient and nec- 
essary conditions for entanglement, we were able to test 
the reliability of our results. The percentage of mislead- 
ing conclusions as a function of the number of samples 
used in the SDP is plotted in Fig. 2. Notice that for 
N > 500 no mistake was made. 




FIG. 2: Percentage of wrong results X number of samples 
(N), for 2x2 (dashed line) and 2x3 (solid line) systems. 

As a third example, we will analyze a three-partite 
bound entangled state derived from the context of 
the micxtendible product bases (UPB) [18]. Con- 
sider the complementary state to the Shifts UPB: 
{|0, l,+),|l,+,0), 1-^,0,1), I-. -,-)}, where ± = (|0) ± 
\1))/V2. We have calculated the OEW for the three bi- 
partite partitions and for the three-partite partition. The 
results of the computation with 2000 samples are sum- 
marized in Table I. Wc can conclude that the state is 
separable with respect to the bipartite splits, whereas 
it is entangled with respect to tripartite product states. 
These same results were obtained using a different ap- 
proach in [18]. 

We have considered so far only the discrimination be- 
tween entangled and separable states. Actually, the 
structure of multipartite quantum entanglement is much 
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Partition 


Ir(WoptP) 


\ yW ) 


\ 


A-BC 


-S.SgarlO"" 


0,063 


-4.34a;10"'' 


B-AC 




0,040 


-S.TSajlO-" 


C-AB 


-1.12a;10-'' 


0,087 


-S.eQajlO-" 


A-B-C 


-3.17a;10~'' 


0,002 


-9.23a;10~' 



TABLE L Results of the method for the three-partite bound 
entangled state complementary to the Shifts UPB 

richer [17]. A n-partite density operator is a in- 

separable state if it is possible to find a decomposition to 
it such that, in each pure state term, at most m parties 
are entangled among each other, but not with any mem- 
ber of the other group of n — to parties. Furthermore, 
even in the class of m-separable states, there exist dif- 
ferent types of entanglement, i.e, states which cannot be 
converted to each other by local operations and classical 
communication protocols (LOCC). Since the subspace of 
m-separable density operators is convex and closed, it is 
also possible to apply the Hahn-Banach theorem to it and 
establish the concept of entanglement witness to (to-|- 1)- 
partite entanglement. In order to do that, consider the 
index set P = {1, 2, n}. Let Si be a subset of P which 
has at most m elements. Then W is an (m -I- l)-partite 
entanglement witness if: 

s,^ (t/-! ® ... ® WW\Tp)s,^ ® ... ® |V')s.„ > 

V S^,,...,S, such that (9) 
ULi^*. =^ and 5.,n^. ={} 

Therefore, it is possible to apply the same methods de- 



veloped earlier to (to -|- l)-partite EW, where one has to 
minimize Tr{Wpi „) subject to the RSDP derived from 

As a final example, we determine a tripartite- 
entanglement OEW for the GHZ state \4>ghz) = 
-^(|000) + I 111)), i.e., an operator which separates 

\ipGHz) fi'om the set of bi-separable density matrices. 
Also in this case, using N — 2000 samples, our proce- 
dure has found a genuine OEW for the state: 

Wopt = i(|001)(001| + |010)(010| + |011)(011| 

-F|100)(100j + |101){101| + |110){110| - |000)(111| 
-|111){000|) 

One can easily check that this matrix is indeed positive 
semidefinite over the separable states. 

In summary, we have constructed a procedure to de- 
termine with arbitrary probability and accuracy optimal 
entanglement witness for every entangled state. Thus, 
considering the NP-hardness of the separability problem, 
this approximate method is of great importance to the 
development of the theory of entanglement. The search 
of others approximate algorithms for the optimization of 
EW with improved performance is an interesting problem 
for further research. 
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